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Abstract 

A consistent interpretation of the Higgs data requires the same precision in the Higgs boson 
masses and in the trilinear Higgs self-couplings, which are related through their common 
origin from the Higgs potential. In this work we provide the two-loop corrections at order 
0{ataa) in the approximation of vanishing external momenta to the trilinear Higgs self¬ 
couplings in the CP-violating Next-to-Minimal Supersymmetric extension of the Standard 
Model (NMSSM). In the top/stop sector two different renormalization schemes have been 
implemented, the OS and the DR scheme. The two-loop corrections to the self-couplings 
are of the order of 10% in the investigated scenarios. The theoretical error, estimated both 
from the variation of the renormalization scale and from the change of the top/stop sector 
renormalization scheme, has been shown to be reduced due to the inclusion of the two-loop 
corrections. 
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1 Introduction 


While the discovery of the Higgs boson by the LHC experiments ATLAS and CMS 
certainly marked a milestone for particle physics, it also triggered a change of paradigm: The 
Higgs particle, formerly target of experimental research, has become a tool in the quest for 
our understanding of nature. Although the Standard Model (SM) of particle physics has been 
tested at the quantum level and the discovered scalar particle behaves SM-like there are 
experimental and theoretical arguments to assume it to be a low-energy effective theory of a 
more fundamental theory appearing at some high scale. In the absence of any direct observation 
of new states the study of the Higgs boson and its properties may reveal the existence of beyond 
the SM (BSM) physics. In particular, the discovered particle could be the SM-like Higgs boson 
of the enlarged Higgs sector of a supersymmetric extension of the SM. Supersymmetric (SUSY) 
theories [4 -18 require the introduction of at least two complex Higgs doublets in order to give 


masses to up- and down-type quarks and ensure an anomaly-free theory. This minimal setup is 
extended by a complex singlet superfield in the Next-to-Minimal Supersymmetric extension of 
the SM (NMSSM) [^-34 . After electroweak symmetry breaking (EWSB) the NMSSM Higgs 


sector features seven Higgs bosons, which in the CP-conserving case are three neutral CP-even, 
two neutral CP-odd and two charged Higgs bosons. In contrast to the Minimal Supersymmetric 
extension (MSSM) [^-38 in the NMSSM CP-violation can occur in the Higgs sector already 


at tree level. The additional sources of CP-violation in SUSY theories are interesting not only 
because they clearly mark physics beyond the SM, but also because CP-violation is an important 
ingredient for successful baryogenesis [^. From a phenomenological point of view it entails a 
plethora of interesting new physics (NP) scenarios not excluded by experiment yet. 

In order to study NP extensions, to properly interpret the experimental data and to be able 
to distinguish different BSM realizations, from the theory side we need as precise predictions 
as possible not only for experimental observable^ but also for the parameters of the theory 
under investigation. In the Higgs sector these are in particular the Higgs boson masses and 
couplings. In the recent years there has been quite some progress in the computation of the 
higher order corrections to the Higgs boson masses of both the CP-conserving and CP-violating 
NMSSM. Thus in the CP-conserving NMSSM after the computation of the leading one-loop 
(s)top and (s)bottom contributions [41| - |45| and the chargino, neutralino as well as scalar one- 
loop contributions at leading logarithmic accuracy 


46 


the full one-loop contributions in the 
and subsequently in 


47 


48 


In 


47 


DR renormalization scheme have first been provided in 
also the order O^atOs + c^bcts) corrections in the approximation of zero external momenta have 
been given, and recently, first corrections beyond order 0{atas -|- aftCts) have been published 
Our group has calculated the full one-loop corrections in the Feynman diagrammatic 

In the CP- 


in 49,50 


51 


approach in a mixed DR-on-shell and in a pure on-shell renormalization scheme 
violating NMSSM the contributions to the mass corrections from the third generation squark 
sector, from the charged particle loops and from gauge boson contributions have been computed 
in the effective potential approach at one loop-level in Refs. (52 -56 . The full one-loop and loga¬ 


rithmically enhanced two-loop effects in the renormalization group approach have subsequently 
been given [^. We have contributed with the calculation of the full one-loop corrections in 


the Feynman diagrammatic approach 58 and recently provided the two-loop corrections to the 


neutral NMSSM Higgs boson masses in the Feynman diagrammatic approach for zero exter¬ 
nal momenta at the order 0{atas) based on a mixed DR-on-shell renormalization scheme 59 


^Neutral NMSSM Higgs production through gluon fusion and bottom-quark annihilation including higher order 


corrections has been discussed in 40 
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Several codes have been published for the evaluation of the NMSSM mass spectrum from a 
user-defined input at a user-defined scale. The Fortran package NMSSMTools [60| - [62| computes 
the masses and decay widths in the CP-conserving Z^-invariant NMSSM and can be interfaced 
with SOFTSUSY |63[[64] , which provides the mass spectrum for a CP-conserving NMSSM, also 
including the possibility of Z^ violation. Recently, it has been extended to include also the 
CP-violating NMSSM 65 . The spectrum of different SUSY models, including the NMSSM, 


can be generated by interfacing SPheno |66, 67 with SARAH 149 68-71 . This is also the case for 


the recently published package FlexibleSUSY (72 73 , when interfaced with SARAH. All these 


codes include the Higgs mass corrections up to two-loop order, obtained in the effective poten¬ 
tial approach. The program package MMSSMCALC (74,75 on the other hand, which calculates 


Higgs self-couplings 79 -81 


the Higgs sector of EWSB 82 - 84 


the NMSSM Higgs masses and decay widths in the CP-conserving and CP-violating NMSSM, 
provides the one-loop corrections and the 0{atas) corrections in the full Feynman diagrammatic 
approach, where the latter are obtained in the approximation of vanishing external momenta. 

The Higgs self-couplings are intimately related to the Higgs boson masses via the Higgs 
potential. For a consistent description therefore not only the Higgs boson masses have to be 
provided at highest possible precision, but also the Higgs self-couplings need to be evaluated 
at the same level of accuracy. The trilinear Higgs self-coupling enters the Higgs-to-Higgs decay 
widths. These can become sizable in NMSSM Higgs sectors with light Higgs states in the 
spectrum (76| - |78( , and via the total width these decays sensitively alter the branching ratios 
of these states. Also Higgs pair production processes are affected by the size of the trilinear 

Their determination marks a further step in our understanding of 
We have provided the one-loop corrections to the trilinear 
Higgs self-couplings for the CP-conserving NMSSM (^. They have been calculated in the 
Feynman diagrammatic approach with non-vanishing external momenta. The renormalization 
scheme that has been applied is a mixture of On-Shell (OS) and DR conditions. In this paper 
we present, in the framework of the CP-violating NMSSM, our computation of the dominant 
two-loop corrections due to top/stop loops to the trilinear Higgs self-couplings of the neutral 
NMSSM Higgs bosons. In addition, we give explicit formulae for the leading one-loop corrections 
at order 0{at). We use the Feynman diagrammatic approach in the approximation of zero 
external momenta and furthermore work in the gaugeless limit. We find that the determination 
of the two-loop corrections reduces the error on the trilinear Higgs self-coupling due to unknown 
higher order corrections and hence contributes to the effort of providing precise predictions for 
NMSSM parameters and hence observables. We have furthermore expanded for this paper the 
full one-loop corrections with full momentum dependence to include CP-violating effects. 

The outline of the paper is as follows. In sectionj^we set our notation, introduce the NMSSM 
Higgs sector and present the determination of the loop-corrected effective trilinear Higgs self¬ 
couplings. Section (^ is then dedicated to the numerical analysis. We discuss the effects of the 
loop corrections on the trilinear Higgs self-couplings and the implications for Higgs-to-Higgs 
decays. Section contains our conclusions. 


2 The effective trilinear Higgs self-couplings in the NMSSM 

In this section we present the details of the calculation of the effective trilinear Higgs self¬ 
couplings at order 0{at) and at order 0{atas). We closely follow the convention and notation 
of our paper on the Higgs mass corrections in the complex NMSSM at order 0{atas) [^. We 
therefore repeat here only the most important definitions relevant for our calculation. We work 
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in the framework of the complex NMSSM with a scale invariant superpotential and a discrete 
symmetry. In terms of the two Higgs doublets and Hu, and the scalar singlet S, the Higgs 
potential reads, 

Vh = (|AS'P + + (|AS'p + rn\j^)HHHu,i + rn\\S\^ 

+ ^-{gl+gl){Hl,Hd, - HHHu,if + ^gl\Hl^Hu,\^ ( 2 . 1 ) 

+1 — XH^^Huj + XAxSHd^Huj + -kA^S^ + h.c.] . 

The indices of the fundamental representation of SU{2)l are denoted by i,j = 1,2, and is 
the totally antisymmetric tensor with ei 2 = = 1. The dimensionless parameters A and k and 

the soft SUSY breaking trilinear couplings A\ and can in general be complex. The U{1)y 
and SU{2)l gauge couplings are given by gi and ( 72 , respectively. In order to obtain the Higgs 
boson masses, trilinear and quartic Higgs self-couplings from the Higgs potential, the Higgs 
doublets and the singlet field are replaced by the expansions about their vacuum expectation 
values (VEVs), Vd,Vu and Vg, 




hi 


Hu = 6**^“ 




ht 

u T hu 


5 


gi^Ps 

--^(us + hg + ittg), (2.2) 


where two additional phases, ipu and (fg, have been introduced. Note, that in order to keep 
the Yukawa coupling neutral we absorb the phase (pu into the left- and right-handed top fields. 


which of course affects all couplings involving only one top quark 59 


We work in the approximation of zero external momenta and call the thus derived self¬ 
couplings ’effective’ self-couplings. The (loop-corrected) self-couplings are automatically real 
in this approach. In the interaction basis, the effective trilinear Higgs self-couplings at order 
0{atag) can be cast into the form 


(2.3) 


with (j) = ihd,hu,hg,ad,au,as) and i,j,k = I,..., 6. The first term represents the tree-level 
trilinear couplings, which can directly be derived from the tree-level Higgs potential Eq. (2.1) 
by taking the derivative 

d^Vn 


Explicit expressions for these couplings can be found in Appendix]^ The second and third terms 
denote the one- and two-loop corrections to the Higgs self-couplings. They can be obtained by 
either taking the derivative of the corresponding loop-corrected effective potential or by using the 
Eeynman diagrammatic approach in the approximation of zero external momenta. At one-loop 
level we use both methods and find that the results obtained in these two different approaches 
agree as expected. However, at two-loop level, for the sake of automatization of our codes we 
solely employ the Feynman diagrammatic approach^ Therefore only the latter is described in 
the following. 


In order to obtain the effective trilinear couplings in the mass eigenstate basis, the self¬ 
couplings in the interaction basis have to be rotated to the mass basis by applying the rotation 


^In the effective potential approach the derivatives which are taken to get the Higgs self-couplings lead to very 
large intermediate expressions, that are not practical to be used for automatization. 
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matrix . In detail, we have, 


(2.5) 


where Z = 1, 2 stands for the loop order and <I> for the loop-corrected mass eigenstates. These are 
denoted by upper case H and ordered by ascending mass with Mh^ < Mh 2 < < Mh_^ < 

Mh^- The neutral Goldstone boson G has been singled out. Note in particular, that the mass 
eigenstates are no CP eigenstates any more since we work in the CP-violating NMSSM. In order 
to be as precise as possible in the computation of the loop-corrected trilinear Higgs self-couplings 
in the mass eigenstate basis, we employ the most precise rotation matrix that is available. 
This means, that we rotate to the mass eigenstates Hi {i = 1,..., 5) at two-loop order. The loop- 
corrected rotation matrix is computed by the Fortran package NMSSMCALC 


74 


75 


where the 

zero momentum approximation is employed, so that the matrix is unitary. In particular, the 
rotation matrix includes the complete electroweak (EW) corrections at one-loop order and the 


order 0{atas) corrections at two-loop level. For more details see [^,58,59 


The rotation matrix can be decomposed in the rotation matrix TZ, that rotates the 
interaction eigenstates to the tree-level mass eigenstates = (/ii,/ 12 , ha,/i 4 ,/is, G), singling 
out the Goldstone boson G, and in the finite wave-function renormalization factor Z, cf. [51[|79|, 


TZ[] — ZijTZjs , 


i,j,s = 1,..,6 . 


( 2 . 6 ) 


With this definition, the loop-corrected effective trilinear couplings between the Higgs bosons in 
the 2-loop mass eigenstate basis are hence given by (i, J, k, k' = 1,..., 5) 




(2.7) 


where the couplings in the tree-level mass eigenstates are obtained from Eq. (2.3) by 

rotation with TZ. 


2.1 The order 0{at) corrections 

In this subsection we present the one-loop corrections at order 0{at). Due to the large top quark 
Yukawa coupling, at one-loop level the corrections from the top/stop sector are the dominant 
corrections to the Higgs boson masses and self-couplings. This is in particular true for the SM- 
like Higgs boson. The latter must be dominantly /i„-like, inducing via the top loop a sufficiently 
large coupling to the gluons, so that its rates are in accordance with the measured signal rates of 
the discovered Higgs boson, which at the LHG is dominantly produced through gluon fusion. The 
restriction to the order 0{at) corrections with large top/stop masses in the loops furthermore 
ensures the approximation of zero external momenta to be reliable. This approximation breaks 
down if the masses of the particles running in the loops are smallj^ For the numerical analysis 
presented in section we took care to choose scenarios where all possibly involved loop particles 
are sufficiently heavy so that not only the approximation of zero external momenta works well 
but also the order 0{a) corrections do not play a significant role. The full EW and the order 
0{at) corrections differ by less than 4% for the chosen scenarios as we explicitly verihed. 

In the following we give the analytic formulae for the one-loop order 0{at) corrections to 
the trilinear Higgs self-couplings in the interaction basis at vanishing external momenta. These 

^Scenarios with light Higgs bosons are mostly precluded as otherwise the kinematically allowed Higgs-to-Higgs 
decays would lower the branching ratios of the SM-like Higgs boson into the other SM particles to values not 
compatible with the experimental data any more. However, other light particles running in the loops could spoil 
the validity of the zero momentum approximation. 
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formulae are compact enough to be easily implemented in computer codes. In order to extract 
only the 0{at) and later on the 0{atas) corrections we neglect all D-term contributions to the 
Higgs potential and to the stop mixing matrix i.e. we work in the gaugeless limit, where the 
electric coupling e and the W and Z boson masses Mw and Mz are taken to be zero but the 
vacuum expectation value v and the weak mixing angle 9w are kept finite. In this approximation, 
the stop mass matrix reads 


Mi = 


m?- + 
Ws 


mt ( — 


tan /3 


tan R 


mi + mi 
tR ^ 


( 2 . 8 ) 


where mt denotes the top quark mass and the effective higgsino mixing parameter 

/^eff — 

and the ratio of the two VEVs Vu and Vd, 


^/2 


tanp = — , 
Vd 


(2.9) 


( 2 . 10 ) 


have been introduced. The soft SUSY breaking masses and m^^ are real, whereas the 
trilinear coupling At = \At\At) is in general complex. The matrix is diagonalized by a 
unitary matrix Ui, rotating the interaction states ti and tR to the mass eigenstates ti and t 2 , 


= Ui {iL,iR)^ 
diag(m|, m|) = Ui uj . 


( 2 . 11 ) 

( 2 . 12 ) 


The order 0{at) corrections to the trilinear Higgs self-couplings in the interaction basis are 
decomposed as 

-b . (2.13) 

The first term denotes the unrenormalized part arising from the one-loop diagrams with tops 
and stops running in the loops. The explicit expressions for are given in Appendix [b| 

The contributions from the parameter counterterms are collected in the second part A^A^!^ ^^. 
Their explicit expressions in terms of the counterterms, defined in the following, are given in 
Appendix [C} 

For the order 0{at) and and the order 0{atas) corrections, we need to renormalize the 
following set of parameters 


59 


0 


thaAhuAhs,taa,tas,Mjj±,V,tan/3, |A| , 


(2.14) 


where tp, (j) = hd,hu,hs,ad,as denote the five independent tadpoles, Mr± stands for the mass 
of the charged Higgs boson and v ~ 246 GeV is given by 


2 2,2 
V =V^ + Vd. 


(2.15) 


^As we work in the gaugeless limit, i.e. e = 0 and Mw = Mz = 0 but v ^ 0 and sin^vc A 0) it is convenient 
to choose V and sin6w in the computation of the higher order corrections, instead of Mw and Mz. Note that 


sin^iv does not appear in the Higgs potential in the gaugeless limit. See also 59 , for more details. 
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In order to renormalize the parameters, they are replaced by the renormalized ones and the 
corresponding counterterms as follows: 


44 with (j) = hd,hu,hs,ad,as , (2.16) 

Mjj± Mjj± + , (2.17) 

V ^ V + 6% + 6% , (2.18) 

tan /3 —)• tan /3 + (l®tan/3 + (l®tan/3 , (2.19) 

|A|^|A|+5ffl|A|+<5^)A| . (2.20) 


Here the superscript (1) denotes the counterterms of 0{at) and the superscript (2) the coun¬ 
terterms of 0{atas). 

In addition to the parameter renormalization, also the wave function renormalization of the 
Higgs helds is needed in order to obtain a UV finite result. At 0{at) and 0{atas), only the 
Higgs doublet Hu has a non-vanishing wave function renormalization counterterm , which is 
introduced as 

^ Hu . (2.21) 

The parameters are renormalized in a mixed OS-DR renormalization scheme as described in 
In this scheme part of the parameters, that are directly related to “physical” quantities, are 
renormalized on-shell, and the remaining parameters are defined via DR conditions, as 

thdi^hui^hsHa^itusi ^• ( 2 . 22 ) 

'-V-' -V-^ 

on-shell scheme DR scheme 


59 


While it is debatable if the tadpole parameters can be called physical quantities, their introduc¬ 
tion is motivated by physical interpretation, so that in slight abuse of the language we call their 
renormalization conditions on-shell. For the wavefunction renormalization of the Higgs fields, 
the DR scheme is employed. Note that this procedure is applied for both the order 0{at) and 
the order 0{atas) corrections. We do not repeat the renormalization conditions here, since they 
are introduced in detail in [59| . We give, however, the explicit expressions for the counterterms 
of order 0{at). For the OS renormalization constants at order 0{at) we find in D = 4 — 2e 
dimensions: 


5% 

V 

+ 

+ 


4 


2s2 


V Ml 


«w 


,y| , 

2 a4 


3?7T'^ 

~ + |4J^Ao(m|J + |44Ao(m|)] 


mt 




I4J + 


Bo(0,mJ^,m?^) 


+ I 1^2 I + 


\X\tpVs\HfJ 

V2 


Bo(0,m|^,m|)) 


3\X\m‘^Vs {V2cfi\X\vs - 2|Af|s^CyJ 4 , ^ 


(2.23) 


(2.24) 

(2.25) 
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with 


and 


= - 


StuP' 1 

* V2\At\\X\vsCpc^^ - 2 si 3 (\At\^ + m? + m? - 2mf 


IhTT^s^^ e L 


3m^ 

IGTr^s^n 


V 2 \At\\X\vsCi 3 C^^Fi -2s^(^\At\‘^Fi + m| + m| - 2rnf 



2 1 

- 

log 


= 

Vs 



SlAtWXlm^VsS 

ip. A , 

8y/2ir'^Sj3V 

(T 


= 



3mt 1 

3 

VS- 

e 

IGvr^n^ 


+\Uij‘^Fo{ml,ml^ 


3mt 1 

3 

Vi 

e 

IGvr^n^ 


Fo{x,y) = 

x + y- 


Fi = 

mf — r 

C2 


m~ 

_^ 

Q2 




2 \ 1 


m. 


rrit - 2mt log ^ \ Fo{m^^, 


- ‘^rnl log ^ , m|) 


2x7/ , X 

- log - , 

x-y y 


— m‘~^ + log log 

2 9 

C2 tl 


And for the DR renormalization constants we get: 

1 




—3mj 


.22-2. > = |tan/35®Z^,„ , 5®|A| = 

QTT^v^ Sin p e 2 2 


(2.26) 

(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 

(2.32) 

(2.33) 

(2.34) 


Here (/>a; = (^^ + + T^Atj with and (pA^ being the complex phase of A, and accord¬ 

ingly of At, has been introduced. Furthermore we use cp = cos/3 etc. The functions AQ{rn?^ 
and Bo(p^,mf, 771 , 2 ) denote the scalar one-point and two-point functions, respectively, in the 
convention of [^, and Q is the renormalization scale. 


2.2 The order 0{atas) corrections 

In order to obtain the order 0{atas) corrections we use the Feynman diagrammatic approach 
in the approximation of zero external momenta. These corrections are composed of 


j. = -I- -I- 


(2.35) 


The first part consists of the contributions from genuine two-loop diagrams. These must contain 
either a gluon or gluino or a four-stop coupling in order to give a contribution of order 0{atas)- 























Some sample diagrams are presented in Fig. ^ In the approximation of zero external momenta 
all two-loop three-point functions can be reduced to the product of two one-loop tadpoles and 


to the two-loop one-point integral which are presented analytically in the literature 86-92 
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Figure 1: Sample of genuine two-loop diagrams contributing to the O(atas) corrections to the trilinear Higgs 
self-couplings between Hi, Hj and Hk {i,j,k = 
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Hk 
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-X «C TT 
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Figure 2: Some representative one-loop diagrams with one-loop counterterm insertion contributing to the 0{atas) 
corrections to the trilinear Higgs self-couplings. 


The second part denotes the contributions arising from the one-loop diagrams 

with top quarks and stops as loop particles and with one insertion of a counterterm of order 
0(as) from the top/stop sector. Some representative diagrams for this set are depicted in Fig.j^ 
The parameters of the top/stop and bottom/sbottom sectors are renormalized at order 0{as). 
The bottom quarks are treated as massless, so that the left- and right-handed sbottom states do 
not mix and only the left-handed sbottom with a mass of contributes. We choose the set of 
independent parameters entering the top/stop and bottom/sbottom sector, that we renormalize, 
to be given by 


mt , rriQ^ , and At . 


(2.36) 


Note that At is in general complex. We renormalize these parameters both in the DR and in the 
OS scheme. The definition of their counterterms can be found in 59]^ According to the SUSY 
Les Houches Accord (SLHA) 


93 


When we choose the OS scheme 


94 


convention, and At are given as DR parameters, 

these parameters need finite shifts for the conversion into OS 
parameters. In the DR scheme on the other hand, the given top pole mass must be translated 


®Note that we work in the CP-violating NMSSM, so that we have trilinear couplings between all hve neutral 
Higgs mass eigenstates. 

®Note that our OS scheme does not take into account terms proportional to e. 
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into a DR mass. These translations according to our conventions are described in detail in 59 

The third part consists of contributions arising from the order 0{atas) counterterms. The 
explicit expressions of in terms of 5%, (5®tan/3, 5®|A| and are 

the same as in the one-loop case after replacing the one-loop by the two-loop counterterms. The 
formulae are given in Appendix For the exact definitions of the two-loop counterterms, we 
refer the reader to 1591. 


Our results have been obtained in two independent calculations. For the generation of the 


amplitudes we have employed FeynArts 95,96| using in one calculation a model file created by 
SARAH (esHTOll^ 


and in the other calculation a model file based on the one presented in [98| 
which has been extended by our group to the case of the NMSSM. The contraction of the 
Dirac matrices was done with FeynCalc 


99 . The reduction to master integrals was performed 


using the program TARCER [100 , which is based on a reduction algorithm developed by Tarasov 
[101, 102 ] and which is included in FeynCalc. We have applied dimensional reduction 103,104 


in the manipulation of the Dirac algebra and in the tensor reduction. In our calculation no 75 
terms appear that require a special treatment in D dimensions, so that we take 75 to be anti¬ 
symmetric with all other Dirac matrices. The cancellation of the single pole and double poles 
has been checked. The results of the two computations are in full agreement. We furthermore 
compared our results in the limit of the real MSSM with Ref. [I05] where the two-loop 0{atas) 
corrections to the MSSM Higgs self-couplings were given, and we found agreement between the 
two computations. 


3 Numerical analysis 
3.1 Scenarios 

For the numerical analysis of the impact of the higher order corrections on the Higgs self¬ 
couplings we made sure to choose scenarios that comply with the experimental constraints. 
In order to find viable scenarios we performed a scan in the NMSSM parameter space. We 
checked the scenarios for their accordance with the LHC Higgs data by using the programs 
HiggsBounds 106-108 and HiggsSignals [109 . The programs require as inputs the effective 


couplings of the Higgs bosons, normalized to the corresponding SM values, as well as the masses, 
the widths and the branching ratios of the Higgs bosons. These have been obtained for the SM 
and NMSSM Higgs bosons from the Fortran code NMSSMCALC j^,75 . A remark is in order for 
the loop-induced Higgs couplings to gluons and photons. The effective NMSSM Higgs boson 
coupling to the gluons normalized to the corresponding coupling of a SM Higgs boson with same 
mass is obtained by taking the ratio of the partial widths for the Higgs decays into gluons in the 
NMSSM and the SM, respectively. The QCD corrections up to next-to-next-to-next-to leading 
order in the limit of heavy quarks 110- 119| and squarks 120, 121j are included. As the EW 
corrections are unknown for the NMSSM Higgs boson decays, they are consistently neglected 
also in the SM decay width. The loop-mediated effective Higgs coupling to the photons has 
been obtained analogously. Here the NLO QCD corrections to quark and squark loops including 
the full mass dependence for the quarks 113,122-127 and squarks [128 are taken into account. 
The EW corrections, which are unknown for the SUSY case, are neglected also in the SM. 

For the numerical analysis of the corrections to the Higgs self-couplings, we have chosen 
two parameter sets that fulfill the above constraints. For both scenarios we use the SM input 
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parameters [129, 130 


a{Mz) = 1/128.962, a^^{Mz) = 0.1184, Mz = 91.1876 GeV, (3.37) 

Mw = 80.385 GeV, mt = 173.5 GeV, mp(mp) = 4.18 GeV . 

In the numerical evaluation, however, we chose to use the running It is obtained by 

converting the , that is evaluated with the SM renormalization group equations at two-loop 
order, to the DR scheme. The light quark masses, which have only a small influence on the loop 
results, have been set to 

niu = 2.5 MeV , = 4.95 MeV , mg = 100 MeV and me = 1.42 GeV . (3.38) 

The remaining parameters differ in the two scenarios. Thus we have: 

Scenario 1: The soft SUSY breaking masses and trilinear couplings are chosen as 


mu^rcR = ^dn,SR = ^Qi .2 = = 3 TeV, = 1909 GeV, 

ttIq^ = 2764 GeV, = 1108 GeV, = 472 GeV, = 1855 GeV, 

\Au,c,t\ = 1283 GeV, \Ad,s,b\ = 1020 GeV, \Ae,f,,r\ = 751 GeV, (3.39) 

|Mi| = 908 GeV, jMal = 237 GeV, \Ms\ = 1966 GeV, 

‘fAa,s,b = = V’Aeect = = ‘fM2 = ‘PM3 = 0 . 

The remaining input parameters are given by 

|A|= 0.374, |k|= 0.162, |^^| = 178 GeV , |/refr| = 184 GeV , 

(^A = = 0 , (fAg^ = TT , tan/3 = 7.52 , Mh± = 1491 GeV . (3.40) 

Scenario 2: For the soft SUSY breaking masses and trilinear couplings we chose 


mun,CR = ^dnrsR = = ^1^,2 = ^^rAr = 3 TeV, = 1170 GeV, 

= 1336 GeV, = 1029 GeV, = 2465 GeV, mf^ = 301 GeV 

\Au,cA = 1824 GeV, \Ad,s,b\ = 1539 GeV, \Ae,f,,r\ = 1503 GeV, (3.41) 

|Mi| = 862.4 GeV, jMsI = 201.5 GeV, \Ms\ = 2285 GeV 

‘fAa,s,b = = TT , = V’Mi = <PM2 = ‘PM 3 = 0 • 

And the remaining input parameters are set as follows. 


|A| = 0.629 , |k| = 0.208 , \A^\ = 179.7 GeV , \fies\ = 173.7 GeV , 

= 0 , (^K = vr, tan/3 = 4.02, = 788 GeV . (3.42) 


We follow the SLHA format, which requires ^eS as input parameter. The values for Vg and (ps can 
then be obtained by using Eq. (2.9). In the SLHA format, the parameters A, k, A^,tan/3 as 
well as the soft SUSY breaking masses and trilinear couplings are understood as DR parameters 
at the scale pR = whereas the charged Higgs mass is an OS parameter. We set the SUSY 
scale Mo to 


• ( 3 - 43 ) 

^For tan /3 this is only true, if it is read in from the block EXTPAR as done in NMSSMCALC. Otherwise it is the 
DR parameter at the scale Mz- 
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OS 

Hi 

H2 

Hs 

Hi 

Hs 

mass tree [GeV] 
main component 

71.14 

hu 

117.49 

hs 

211.12 

(Ig 

1491 

a 

1492 

hd 

mass one-loop [GeV] 
main component 

98.65 

hg 

139.17 

hu 

217.27 

(Xg 

1490 

a 

1491 

hd 

mass two-loop [GeV] 
main component 

94.68 

hg 

125.06 

hu 

217.32 

dg 

1490 

a 

1491 

hd 

DR 

Hi 

H2 

Hs 

R4 

Hs 

mass tree [GeV] 
main component 

71.14 

hu 

117.49 

hg 

211.12 

dg 

1491 

a 

1492 

hd 

mass one-loop [GeV] 
main component 

91.60 

hg 

120.00 

hu 

217.36 

dg 

1491 

a 

1491 

hd 

mass two-loop [GeV] 
main component 

94.41 

hg 

124.24 

hu 

217.33 

dg 

1490 

d 

1491 

hd 


Table 1: Scenario 1: Masses and main components of the neutral Higgs bosons at tree and one-loop level and at 
order Oiottcts) as obtained by using OS (upper) and DR (lower) renormalization in the top/stop sector. 


The resulting supersymmetric particle spectrum from the thus chosen parameter values is in 
accordance with present LHC searches for SUSY particles [131 -145 . Note, that in the following 


we will drop the subscript ’eff’ for Furthermore, whenever we will use the expressions OS and 
DR these refer to the renormalization in the top/stop sector. 


3.2 Results for the loop-corrected self-couplings 

In this and the following subsection we discuss the impact of the order 0{atas) corrections on 
the trilinear Higgs self-couplings and on Higgs-to-Higgs decay widths. We start by discussing 
the results for the parameter set called scenario 1 in the previous subsection. The masses of the 
Higgs bosons and their main composition in terms of singlet/doublet and scalar/pseudoscalar 
components at tree level, one-loop and two-loop order are summarized in Table[^both for the OS 
and the DR renormalization in the top/stop sector. The tree-level stop masses in this scenario 
are rather heavy and given by 

OS : mr = 1992 GeV , mr = 2820 GeV , 

- H ’ t2 ’ (3.44) 

DR : = 1911 GeV , = 2768 GeV . ^ ’ 

and for the DR top mass we have = 136.34 GeV. For definiteness, with respect to the mass 
corrections one-loop means here and in the following that we include the full EW corrections 
at non-vanishing external momenta, while at two-loop level the order 0{atas) corrections are 
computed at vanishing external momenta. As can be inferred from the table, the masses of the 
three lightest scalars are substantially different, so that mixing effects due to GP-violation for 
non-vanishing phases cannot be expected to be significant. The reason for choosing this scenario 
are higher order corrections to the trilinear Higgs self-coupling of the SM-like Higgs boson which 
are rather important for this parameter point. This boson is given by the state with the largest 
hu component and a mass value around 125 GeVAt tree level it is the lightest Higgs boson 

® A rather large hu component is required in order to reproduce the experimentally measured production rates. 
They are mainly due to gluon fusion, which is dominantly mediated by top loops for small values of tan/3. 
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Figure 3: Scenario 1: Upper panels: Trilinear self-coupling of the hu interaction state (left) and the mass 
eigenstate h (right) as a function of including the one-loop correction (blue/two outer lines) and also the 
two-loop corrections (red/two middle lines). In the top/stop sector either the OS scheme (solid lines) or the DR 
scheme (dashed lines) has been applied. Lower Panels: Absolute value of the relative deviation of the correction 
using OS or DR renormalization in the top/stop sector, A = |~ hu,h), in percent 
as a function of at one-loop (blue/upper) and two-loop (red/lower). 


Hi that is mainly h^j-like, and its mass thus receives large corrections which are dominantly 
stemming from the top/stop sector. The large corrections shift the Hi mass above the one 
of H 2 so that the two Higgs bosons interchange their roles, as they are ordered by ascending 
mass. At one- and two-loop level it is therefore the second lightest Higgs boson, which is h^j-like. 
For convenience, we denote in the following the mass eigenstate that is dominantly h^-like, by 
h. Furthermore, when we perform comparisons in the interaction basis at different loop-levels, 
these will be done for the hu state. 


In Fig. we show the dependence of the one- and two-loop corrections to the Higgs self¬ 
coupling on the DR parameter At in the two different renormalization schemes applied in the 
stop sector. The one-loop corrections have been obtained at order 0{at) for vanishing exter¬ 
nal momenta. We explicitly verified, that the differences between the one-loop result in this 
approximation and the one including the full one-loop corrections for non-vanishing momenta 
at the thresholc0 are below 4% for the investigated parameter points. Two-loop corrections 
always refer to the order 0{atas) corrections at vanishing external momenta. The left plot of 
Fig.shows the corrections to the self-coupling of hu in the interaction basis, Xhuh^hu- Figure 
(right) displays the loop-corrected self-couplings after rotation to the mass eigenstate h with 
dominant hu component. The rotation to the mass eigenstates is performed with the mixing 
matrix defined in Eq. (2.5) for both the one- and the two-loop curves in the plot. The mass 
values and mixing matrix elements have been computed with NMSSMCALC. Note that at two-loop 
order the hu dominated state is given by the second lightest Higgs boson H 2 , cf. Table The 
dependence on At is more pronounced after rotation to the mass eigenstates. Overall, however, 
the size and shape of the corrections both in the interaction and in the mass eigenstates are 
comparable. At the parameter point of scenario 1 the tree-level coupling Xh^huku = 101.70 GeV 


®The non-vanishing momenta at the threshold have been set to pi = pi = “m-h for two of the external momenta 
and to Pi = 4m| for the remaining one. Here denotes the two-loop corrected mass value of the SM-like Higgs 
boson. 
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Figure 4: Scenario 1: Same as 



but now as a function of • 


in both renormalization schemes. In the OS scheme the one-loop correction increases it by 140% 
while it is decreased by 24% to two-loop order. In the DR scheme the increase is of 74% going 
from tree- to one-loop order supplemented by another increase of 9% when adding the two-loop 
corrections. The reason, why the one- and two-loop corrections differ much more in the OS 
scheme than in the DR scheme can be understood as follows. In the DR scheme the top quark 
mass, which according to the SLHA accord is an OS parameter, has to be converted to the DR 
value. Thereby, the finite counterterm to the top mass, which in the OS scheme is included at 
two-loop level, is already induced at one-loop level in the value of the DR mass. In this way 
some corrections of order 0{atas), which in the OS scheme only appear at the two-loop level, 
are moved to the one-loop level, cf. also 


59 


The lower panels of Fig. display the difference in the self-couplings when using the two 
different renormalization schemes in the top/stop sector. 


A = ' , (3.45) 

where H both refers to the hu dominated mass eigenstate h, and to the hu interaction eigenstate. 
This value gives a rough estimate of the theoretical error in the Higgs self-coupling due to the 
unknown higher order corrections. In the interaction eigenstate it is of order 0(50%) at one-loop 
level, decreasing to roughly 4% at two-loop level. In the mass eigenstate it is about 5% higher 
at both loop orders. The inclusion of the two-loop corrections hence substantially decreases the 
theoretical uncertainty. 

Figure shows the same as Fig. [^but now as a function of the phase ifAt- All other CP- 
violating phases have been kept to zero. The figure shows that the dependence of the loop 
corrections on the phase is almost negligible, as expected for radiatively induced CP-violation. 
The size of the loop corrections and the remaining theoretical uncertainty are of the same order 
as for the variation of At. 

We now turn to the discussion of scenario 2. The masses and dominant composition of the 
mass eigenstates at tree level, one- and two-loop order are summarized in Table In the OS 
scheme again the composition of the mass ordered states changes when going from tree level to 
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OS 

Hi 

H 2 

Hs 

Hi 

Hs 

mass tree [GeV] 
main component 

79.15 

hs 

103.55 

hu 

146.78 

(Is 

796.62 

hd 

803.86 

a 

mass one-loop [GeV] 
main component 

103.45 

hs 

129.15 

Qs 

139.83 

hu 

796.53 

hd 

802.94 

a 

mass two-loop [GeV] 
main component 

102.99 

hs 

126.09 

hu 

128.94 

ds 

796.45 

hd 

803.07 

a 

DR 

Hi 

H 2 

Hs 

Hi 

Hs 

mass tree [GeV] 
main component 

79.15 

hs 

103.55 

hu 

146.78 

ds 

796.62 

hd 

803.86 

a 

mass one-loop [GeV] 
main component 

102.80 

hs 

120.52 

hu 

128.80 

ds 

796.36 

hd 

803.09 

a 

mass two-loop [GeV] 
main component 

103.09 

hs 

124.55 

hu 

128.91 

ds 

796.36 

hd 

803.03 

a 


Table 2: Scenario 2: Masses and main components of the neutral Higgs bosons at tree and one-loop level and at 
order 0(a.taa) as obtained by using OS (upper) and DR (lower) renormalization in the top/stop sector. 


one-loop level and from one- to two-loop level. In contrast to scenario 1 the masses of H 2 and 
H'i are now much closer together, in particular after inclusion of the two-loop corrections. We 
therefore expect CP-violating effects to be more important here. The H 2 state is identified with 
the discovered Higgs boson. The stop masses are again rather heavy with 

OS : mr = 1145 GeV , m,- = 1421 GeV , 

- H ’ t2 ’ (3.46) 

DR : = 1126 GeV , = 1387 GeV . ^ ’ 

In Fig. [^we show the dependence of the Higgs self-coupling of the hu state in the interaction 
basis (left) and of the /lu-like mass eigenstate h (right) at one- (dashed) and two-loop order 
(full) for DR renormalization in the top/stop sector as a function of the phases ^m^-, ^At and 
(Pfj,. For illustrative purposes we have varied the phases in rather large ranges although they 
might already be excluded by experiment. We start from our original GP-conserving scenario 
and turn on the phases one by one. Note, that has been varied such that the CP-violating 
phase <Py = (Pk — ^\ + ‘^P>s — P>u-, that appears already at tree level in the Higgs sector, remains 
zero, ie. ^px and ips were varied at the same time as p)\ = 2p>s = 2 / 3 (^^, while (/?«, and p>u are 
kept zero. As expected, the loop-corrected couplings show a somewhat larger dependence on 
pi At than in scenario 1, in particular in the mass eigenstate basis. Dehning as 

"""" = —wm— 

we have in the mass eigenstate basis H = h the variations 

= <a = 1.6% and = 2.7% (3.48) 

for the two-loop corrected self-coupling. Note, that the one-loop corrected self-couplings show 
a dependence on the phase of M 3 , although the genuine diagrammatic gluino corrections only 
appear at two-loop level. This dependence enters through the conversion of the OS top quark 


(3.47) 
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Figure 5: Scenario 2: Upper panels: Trilinear self-coupling of the hu interaction state (left) and the mass eigenstate 
h (right) at one-loop (dashed lines) and two-loop (solid lines) as a function of the phases Pfi, (green/grey), ipAt 
(red/black upper) and Pm^. (blue/black lower). Lower Panels: Size of the relative correction of order to the 
Higgs self-coupling with respect to the (n— 1)‘^ order - i.e. A = ~ percent for H = hu 

(left) and H = h (right) as a function of the phases (green/grey), pAt (red/black) and Pm^ (blue/black) for 
n = 2 (solid line) and n = 1 (dashed line). The red and blue lines almost lie on top of each other. In the top/stop 
sector we have applied DR renormalization. 


mass to the DR mass. Overall, the dependence of the loop corrected self-couplings on the 
CP-violating phases is smaller in the interaction states than in the mass eigenstates, which are 
obtained by rotating the interaction states with the mixing elements obtained from the loop 
corrected masses, that also depend on the CP-violating phases. 

The lower panels show the relative corrections, defined at order n = 1,2 as 

|,(n) _ 

A = ^hhI (349^ 

In the interaction basis they are of order ~ 70 — 80% for the one-loop corrections relative to the 
tree-level coupling and are somewhat larger than the corresponding values in the mass eigenstate 
basis, which are of order ~ 50 — 60%. For the two-loop coupling relative to the one-loop coupling 
the corrections are signihcantly reduced to about 5 — 8% in both the interaction and the mass 
eigenstate basis. The two-loop corrections hence considerably reduce the theoretical uncertainty. 

In order to further study the theoretical uncertainty, we show in Fig. [^for scenario 1 (left) 
and scenario 2 (right) the scale variation of the trilinear Higgs self-coupling in the mass eigenstate 
h at one- and at two-loop order. We have varied the renormalization scale between 1/3 and 
3 times the central scale /tq = Mg. The scale variation affects the DR parameters entering the 
calculation. In the absence of an implementation of the 2-loop renormalization group equations 
(RGE) for the complex NMSSM, which is devoted to future work, we obtain the parameters at 
the different scales by exploiting the relation between DR and OS parameters, as explained in 
Appendix [P) This should approximate the results obtained from the RGE running rather well, 
in case the scale is not varied in a too large range. Since the scale variation provides only a rough 
estimate of the error made by neglecting higher order corrections this approach is sufficient for 
our purpose. As can be inferred from the figures, in scenario 1 the one-loop coupling is altered 
by up to 7% compared to its value at the central scale in the investigated range. This reduces 
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Figure 6: Scale dependence of the trilinear coupling Xhhh at one- (blue/dashed) and two-loop order (red/full) 
for scenario 1 (left) and scenario 2 (right). The scale fiu = ^Ms has been varied in an interval of ^ = 1/3...3 
around the central scale /xq = M^. The lower plots show the variation in percent compared to the central scale, 
i-e. A = [XhhhiiJ'R) — XhhhilO'O)]/XhhhilJ'O) ■ 


to 2-5% at two-loop order. In scenario 2 the corresponding numbers at one-loop order are 3.5% 
compared to up to 2.5% for the two-loop coupling. As expected, the scale dependence reduces 
when going from one- to two-loop order. Note, however, that these numbers should not be taken 
as estimate for the residual theoretical uncertainty. 


3.3 Phenomenological implications 

We now turn to the discussion of the phenomenological implications due to the loop-corrected 
Higgs self-couplings. Higgs self-couplings are involved in Higgs-to-Higgs decays and in Higgs 
pair production processes. At the LHC, pair production dominantly proceeds through gluon 
fusion. This process, however, includes EW corrections beyond those approximated by the loop- 
corrected effective trilinear couplings. As they are not available at present we will not discuss 
Higgs pair production further and concentrate on Higgs-to-Higgs decays. 

The decay width for the Higgs-to-Higgs decay Hi —)• HjHk including the two-loop corrections 
to the Higgs self-coupling is obtained from 


V{H, ^ HjHk) 


167r/ Mfj^ 




(3.50) 


where / = 2 for identical final state particles and / = 1 otherwise. The decay amplitude is 
denoted by and A = {x — y — — ^yz is the two-body phase space function. In case 

of CP-violation all Higgs-to-Higgs decays between the five neutral Higgs bosons are possible, 
if kinematically allowed, so that i,j. A: = 1, ...,5. In the CP-conserving case, however, only the 
trilinear Higgs couplings between three CP-even or one CP-even and two CP-odd Higgs bosons 
are non-vanishing. The matrix element is given by 


5 


(3.51) 
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Here is the loop corrected trilinear Higgs self-coupling in the tree-level mass eigenstate 

basis, where the Goldstone boson has been singled out. We here include at one-loop level the 


full electroweak corrections at 7 ^ 0, where we set the 4-momenta of the external Higgs 
particles equal to the respective loop-corrected Higgs mass values ^ = M^, , ^ as obtained 
with NMSSMCALC and at two-loop the order O(atas) corrections at p^ = 0 The 

proper on-shell conditions of the external Higgs bosons as required in the decay process are 
ensured by rotating the tree-level mass eigenstates to the loop corrected mass eigenstates 

^i,j,k with the matrix Z, cf. 51,79 . In this calculation we include at one-loop order the 
full electroweak corrections at non-vanishing external momenta. At two-loop order as usual 
the order 0{atas) corrections which are available only at = 0 are taken into accountThe 
accounts for the contributions stemming from the mixing of the CP-odd components 
of the external Higgs bosons with the Goldstone and with the Z boson, respectively. These 
contributions, which are evaluated by setting the external momenta to the tree-level masses 
in order to maintain gauge invariance, are small already at one-loop order compared to the 
remaining contributions to the decay amplitude, as has been shown in [^. We hence do not 
include the two-loop contributions, that can safely be expected to be negligible. 

In the plots below we show apart from the two-loop corrected decay widths also the ones at 


one-loop order. The only change required to adapt formula (3.51) to this case is in Th.,h-,hi., 
where solely the one-loop corrections to the vertex functions together with the corresponding 
counterterms are included. In particular we also use the two-loop corrected mass eigenstates 
and mixing matrix elements for the external particles (apart from the mixing contribution with 
the Goldstone and Z boson of course). 

The scenario which the following discussion is based on and which has been checked to be 
compatible with the constraints from the LHG Higgs data, is given by: 

Scenario 3: The soft SUSY breaking masses and trilinear couplings are chosen as 


= 3 TeV, = 1940 GeV, 

= 2480 GeV, = 1979 GeV, = 2667 GeV, ruf^ = 1689 GeV, 

= 1192 GeV, \Ad,s,b\ = 685 GeV, \Ae,R,r\ = 778 GeV, (3.52) 

|Mi| = 517 GeV, IM 2 I = 239 GeV, IM 3 I = 1544 GeV, 

‘PAa,s,b = = 0 > ^Au,c,t = ^ = <PM2 = = 0 . 

The remaining input parameters are given by 


|A|= 0.267, |«:|= 0.539, |A| = 810 GeV , |peff| = 104 GeV , 

= V’k = =‘Pu = 0 , (pA^=Tr, tan/3 = 8.97, = 613 GeV . (3.53) 


This results in the Higgs mass spectrum given in Tableat tree, one- and two-loop level together 
with the main singlet/doublet and scalar/pseudoscalar components at each loop level. 

In Fig. (left) we show the partial decay width for the decay of the heavy h^-like Higgs 
boson H 4 into a pair of SM-like Higgs bosons H 2 , including the higher order corrections to the 
Higgs-self-couplings at one- and two-loop level as obtained from Eq. (3.50). We start from the 


the loops the tree-level masses for the Higgs bosons are used to ensure the cancellation of the UV diver¬ 
gences. 

^^As we investigate here the decay of heavy particles in the initial and final states, it makes sense not to work 
in the zero momentum approximation if possible and include at one-loop level the full momentum dependence. 
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OS 

Hi 

H 2 

Hs 

Hi 

H 5 

mass tree [GeV] 
main component 

49.17 

hs 

99.83 

hu 

609.21 

a 

611.77 

hd 

715.92 

(Is 

mass one-loop [GeV] 
main component 

87.36 

hs 

139.10 

hu 

608.71 

0 

611.37 

hd 

694.73 

ds 

mass two-loop [GeV] 
main component 

83.66 

hs 

124.95 

hu 

608.73 

0 

611.37 

hd 

694.76 

ds 

DR 

Hi 

H 2 

Hs 

Hi 

H 5 

mass tree [GeV] 
main component 

49.17 

hs 

99.83 

hu 

609.21 

a 

611.77 

hd 

715.92 

ds 

mass one-loop [GeV] 
main component 

80.66 

hs 

119.68 

hu 

608.72 

0 

611.37 

hd 

694.79 

ds 

mass two-loop [GeV] 
main component 

83.03 

hs 

124.34 

hu 

608.71 

a 

611.36 

hd 

694.78 

ds 


Table 3: Scenario 3: Masses and main components of the neutral Higgs bosons at tree and one-loop level and at 
order Oiottas) as obtained by using OS (upper) and DR (lower) renormalization in the top/stop sector. 


parameter point of scenario 3 with vanishing phase (/^^ = 0. The phase is then varied in the 
range —7r...7r, such that at tree level the CP-violating phase (py in the Higgs sector vanishes. As 
expected the dependence of the decay width on the CP-violating phase induced through the loop 
corrections is small, remaining below the per-cent level. For = 0 the tree-level decay width 
in the OS scheme is 0.171 GeV and 0.186 GeV in the DR scheme)^ In the latter the one-loop 
corrections increase the decay width by 6.5% and the two-loop corrections add another 2.0% on 

^^Note, that of course also in the tree-level decay width we use the Hu and H 2 mass values including the 
two-loop corrections and rotate to the mass eigenstates with the corresponding mixing matrix elements. 



^ 4 . 8 - '' '' 

4.6 |- ... - 

-TT -7r/2 0 4-7r/2 -Htt 

‘Pf. 





Figure 7: Scenario 3: Upper: Loop-corrected decay width (left) and branching ratio (right) of the Higgs-to-Higgs 
decay H 4 —> H 2 H 2 as a function of with the top/stop sector renormalized in the OS (blue/two outer lines) and 
in the DR scheme (red/two inner lines) at 1-loop (dashed) and at 2-loop (full) order. Lower: Relative deviation 
between the two renormalization schemes A = |0(DR) — 0(0S)|/0(DR) for O = F (left) and O — BR (right) at 
1-loop (dashed) and 2-loop (full). 
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top of that. In the OS we find a 21% increase at one-loop and a 7% decrease at two-loop so that at 
two-loop order the results in the two renormalization schemes approach each other, as can also be 
inferred from the lower left plot, which shows that the dependence on the renormalization scheme 
decreases from one- to two-loop level. Depending on the renormalization scheme different input 
parameters of the top/stop sector have to be converted to match the required renormalization 
scheme. The corresponding induced two-loop corrections into the one-loop corrections lead to 
a different dependence on the CP-violating phase of the two schemes at one-loop level. At two- 
loop level this difference in the phase dependence is then almost washed out and the scheme 
dependence is about 4.58% independent of For the computation of the branching ratio of the 
decay, shown in Fig. (right), we replace in the program package NMSSMCALC the tree-level decay 
widths T(Hi —7- HjHk) with our loop-corrected ones. The branching ratio, which with 0(3%) is 
very small shows the same trend as the decay width with respect to the loop corrections. 

The non-vanishing CP-violating phase induces through the higher order corrections CP mix¬ 
ing in the Higgs mass eigenstates, such that otherwise not allowed decays of e.g. the CP-odd 
doublet-like into a pair of SM-like H 2 bosons are possible. The branching ratio remains, 
however, tiny, reaching at most 0.58 per mille for ~ 7r/2 in our scenario. 

4 Conclusions 

The search for New Physics and the proper interpretation of the experimental data requires from 
the theoretical side precise predictions of parameters and observables. In this work we computed 
the two-loop corrections to the trilinear Higgs self-couplings in the CP-violating NMSSM. Orig¬ 
inating from the Higgs potential the Higgs boson masses and self-couplings are related to each 
other. For a consistent interpretation therefore the level of accuracy of the self-couplings has to 
match the one of the masses, that have been provided previously up to the two-loop level. Here, 
the two-loop corrections to the self-couplings have been calculated at the same order 0 {atas) at 
vanishing external momenta. We have allowed for two renormalization schemes in the top/stop 
sector, namely OS and DR renormalization. Depending on the scenario and the renormalization 
scheme, the two-loop corrections are of the order of 5-10% relative to the one-loop couplings, 
compared to up to 80% for the one-loop corrections relative to the tree-level values. The investi¬ 
gation of the remaining theoretical uncertainty performed by varying the renormalization scheme 
of the top/stop sector or by changing the renormalization scale confirmed that the theoretical 
error is reduced through the inclusion of the two-loop corrections. As expected the dependence 
on the CP-violating phase due to radiatively induced CP-violation is small and of the order of 
a few percent. 

The trilinear self-couplings are relevant in Higgs pair production processes and Higgs-to-Higgs 
decays, which now become accessible at run 2 of the LHC. While Higgs pair production requires 
the inclusion of further higher order corrections beyond the loop-corrected Higgs self-couplings 
provided in this work, the inclusion of the radiatively corrected trilinear Higgs self-coupling 
improves the prediction for the Higgs-to-Higgs decay rates and related branching ratios. For the 
investigated scenario and decay we find that the two-loop corrections alter the decay width by 
2% (7%) in the DR (OS) scheme with respect to the one-loop level, which is to be compared to 
about 7% (21%) when going from tree- to one-loop level. The dependence on the renormalization 
scheme is reduced from ~ 4.8 — 5.4% in the investigated range of the phase at one-loop level 
to ~ 4.5% at two-loop level. The behaviour in the branching ratio is similar to the one of the 
decay width. 
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In summary, the inclusion of the two-loop corrections at 0{atas) in the approximation of 
vanishing external momenta in the trilinear Higgs self-couplings of the CP-violating NMSSM 
Higgs sector is necessary to match the available precision in the Higgs masses and to allow for 
a consistent interpretation of the Higgs data. Being of the order of 10% they have been shown 
to further reduce the theoretical error due to missing higher order corrections. 
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Appendix 

A Tree-level trilinear Higgs self-couplings 


In this appendix we present the tree-level trilinear Higgs self-couplings in the interaction eigen¬ 
states, Xijk, i,j,k = 1, ...,6 with the correspondences 1 = hd, 2 = hu, 3 = /i*, 4 = o^, 5 = a^, 
6 = tts- They are symmetric in the three indices. Using the short-hand notations Cx = cosx etc. 
and 


i'x — + ^Ax : 

t>y — 2(^5 -|- ip\ ifi^ , 

Pz = 3 (Ps + ‘PAk + ) 


(A.54) 


we have 


Am — 




An2 — — 


I l\|2 


+ |A| SjdV , Alls — |A| Vs, Aii 4 —0, Ans —0, 


Ai16 — 0 , Ai22 — — 


, im2. 1 _ I^a||A|c<^^ 


+ \X\Ci3V, Ai23 = -- 


v/2 


- \X\\K\VsC^y , Ai24 = 0, 


Ai 25 — 0 ) Ai 26 — —-|A||k|usS0j^ , Aiss — |Apuc/3 — |k||A|us^C0j^ , A134 —0, 


Ai35 = x|A||k|UsS 0 j, , Ai36 = —|A||k|uS(/,j^S/3 , Ai44 = Ai45 = 0, 


Ai 46 — 0 ' Ais.'S — — 
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VpMz I Im2^ 1 _ I^a||A|c^^ 


+ |A| CfjV , Ai56 — 
SM'^sis 


V2 


- A KV, 


^‘Py 


A 166 — |A| -I-|«:||A|us/ 3 C(^j^ , A 222 —-, A 223 — |A| Vg, A 224 —0, A 225 —0, 


A226 — 0 , A233 — —|A||k|uc/ 3 C(^j^-I- |A| Sj 3 V , A234 — 2 ’ A235 — 0 , 

A236 = —, A244 =-+ l-^P'^'5/3 > A245 = 0, 

X |2tA||A|c0^ M'^sp 

A246 = -7=- \X\\k\VsC4, , A255 = - , A256 = 0, 

V2 V 

A266 = \M\v^\vci 3 C^y |Aps/3U, A333 = 6 |k|^Us -h V 2 \A^\\k\c^^ , A334 = |k||A|i;s/3S</,^ 
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^ I MM ^ i\\MsfiCpS^y 2 

A335 = \i<.\\X\vci 3 S^y , A336 = -, A344 = |A| Vs, 

Vs 

A345 = , A346 = —\Hi\\X\vSj3C^y , A355 = |A|^t>s , 

X 356 = —\l^\\X\vCj 3 Cfj,y,X 3 QQ = —V 2 \Afi\\K\C(i)^+ 2 \K\'^Vs, A444 = 0 , A445 = 0 , 

3 

A446 0 , A455 0 , A45g 2 11 X\vsS(py , A4()g , 

, ,, , 3 |fi;||A|cflSflSfli„f^ 

A555 = 0, A556 = 0, A566 = —\>i\\X\vci3S^y, Aeee =-■ 


(A.55) 


B The order 0{at) corrections to the trilinear Higgs self-couplings 


The order 0{at) one-loop corrections introduced in Eq. (2.13), {i,j,k = 

1, 6 ), with the same correspondences as introduced in Appendix]^ can be cast into the form 

= -‘^Cpnityl [Fi, mrh^jhi + h\{h^^*h\ + h\h]{h\r) + + c.c.] (B.56) 

- Cpyl [-Fsxyhihhyhjhhyh^hii + F2xyha2i2yhji2iiyhkhi2 + Permutation[z, j, k]^ 
yhii\i\yhji\i\yhtci\i\ 


m 




m 


t2 


- Cpyi 


m 


Fix {yhki2iiyhihjiii2 F /02 




log -^yHt2t2yh,hit2i2 + {k^i) + {k^ j) 


where 


Cf = 


Ihvr^ ’ 


yt = 


V^mt 

vsp 


mr 


, Fix — 2 log -|- 1 , F 2 x — 


mi — mi — mi log 

tl *2 q 

-(B-57) 

[mi — mi Y 

^ tl t2^ 


Fsx — 


mi — mi — mi log 

tl 'V'-' 


h- 

*1 


(m~ — m~ 

tl t 2 


; — 


9 9 

Tn,~ 7n~ 

*2 ’"'t 2 m| - m| - log + m| log 


(mi — mi ) 

n *2 


and the non-vanishing couplings read 
1 . + i 


yh2int 


y hFnin 




v /2 


-h 


v /2 


+ V2mtUi^^Ul^ + ^J2mtUl^^U 


(B.58) 


(B.59) 
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2 2 
yhiinim ~ ~ ’ 
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y h^ifiim 

y hQtnim 
yh 1 /13 in im 

y h 1 Hq in im 
y h2h,2inim 
y h^h^inim 

y h^hQinim 
y h^h^inim 




^/2 


V 2 


2 2 

--X*e-^'^‘Uf W - -Xe^^^Ur W , 

Ux*e-^^‘Ur m - Uxe^^^Ur TM , 
^zA*e-*^«Z^r W - hxe^^^Ur W , 

s*".™,"!.. + y‘Ki...uL ■ 


C The trilinear Higgs self-coupling counterterms 


Here we summarize the one- and two-loop (/ = 1,2) non-vanishing counterterms that 

arise in the computation of the loop-corrected Higgs self-couplings Xijk {i,j,k = 1,...,6). They 
are given in the interaction basis, and read in terms of the various tadpole, mass, wave function 
and parameter counterterms as 


A«AfT 

A«AfT 

A^A^l 


A«AfT 




aWa^I 

A«Af3^e 

A«AfT 

A«AfT 


= 2us/ 3|A|(5^^^|A| -F uc^lApJ^^Han/? -h si3\X\‘^S^% + ^vsi3\X\‘^S^^^Zh^ , 

= 2vs\X\S^X\, 

= 2uc/ 3|A|5('^|A| - uc|s_a|Ap(5('Han^ -F Ci3\X\‘^5^% + vCj3\X\'^6‘'^^Zh^, 

_( I II ^^C;3S/3|A|^^ ^ f ;^ i^i C2/3C^(2M|^±-bu2|A|2),5Wtan/3 

“V 2^'*'^^-'''' Vs ) ' ' 2U. 

CjsSjsS’y'^Mjj^ vciiSii\X\^5^% 

VVs Vs Vs 

_ (s 2 / 3 M|± +vlc^y\n\\X\ +v‘^cpsp\X\^)5^^^Zh^ 

AVs 

= - \vsS^ \K\\X\5^^'^Zh^ -b 

2 " 4 vVsSji 

= v{-c^ySi3\K\ -b 2 c/3|A|)5('^|A| - vc^p\X\{ci3C^y\K\ + S/3|A|)(5(*Han/3 

+ \>'\{-C 4 ,ySp\K\ + c^|A|)5(')u, 

= Ivs\k\s^J^X\ + jVs\K\\X\s^j‘'^'^Zh^ + 

^ 4r 'V'Vg S 

= -vsi3S^Jk\6 ^^]X\ - uc^s^jK||A|5(^Han/? - Si3S^jK\\X\d^%, 

= aWa^t,, AWA^^e = -A«Afi, 

= v{c^ySi3\K\ + 2c/ 3|A|)(5(')|A| -b uc^IAKc/^c^^^IkI - S0\X\)6^^han(3 

+ \Mic<PySi3\K\ -bc/3|A|)5('^u, 


(C.60) 
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= 2 vs\X\ 6 ^^]X\+Vs\X\^ 6 ^^%^, 

= {-vci3C^Jk\ + 2us^|A|)(5(')|A| + vcp\X\{c^ySi3\K\ + C;3|A|)(5('Han/3 
+ \X\{-C 0 C^Jk\ + S/3|A|)<5('V + ^v\X\{ -C 0 C^Jk\ + S 0 \X\) 6 ^^^Zh^, 

A^a^t = a(OaCT, 

= -vc0\k\s^J‘'^]X\ + vcls0\K\\X\s^J^^han/3 - C0 \k\\X\s^^ 6% - ^vc0\K\\X\s0j^^'>Zh^, 

A«aCT = A«AfT, 


A(')AC/e = 


■2^sC4>^\k\ 


C/3g/3|A| \ (;v C2^c|( 2M|± +^^|Ap)jWtan/3 

Vs y ' ' 2 vs 


^ VC 0 S 0 \X\H(-% ^ 00806 ^ 1 ^ 


VVs 


+ 

aWa^t = 
+ 

A«AC3^3 = 
A^A^t = 
A^'^A^l = 
A^'U^^e = 


{s20Mjj± - 2 ,v‘Ic^^\k\\X\ + v‘^C0S0\X\'^)5'^^'>Zh^ 

Avs 

vic 0 C^y\K\ + 2s^|A|)5(^^|A| + vc‘ 0 \X\{-c^yS 0 \K\ + C 0 \X\) 6 ^’-han /3 

\M{c 0 C^y\K\ + s^|A|)(5(% + ^v\X\{c 0 C^y\K\ + S 0 \X\)S^^'^Zh^, 

3v‘^C0S0S0,yt^jK\64X\ 3v‘^4c20S0yt0jK\\X\5^’‘han/3 ^ 6vc0S0S^yt^jK\\X\6^% 

Vs Vs Vs 

i;s/3S0Jk|( 5(')|A| + ?;c^s^Jfi;||A|(^(^Han/3 + S 0 S^Jk\\X\ 6 ‘^%, 

vc0S<0Jk\54x\ - vc0S0S^jK\\X\6‘^’-hanP + C0 S^Jk\\X\ 6^% + ^vc0S^jK\\X\6^'''^Zh^, 
3v‘^C0S0S0y\K\64X\ 3v‘^c^^C20S^jK\\X\5‘'''han P 6 vc0S0S^Jk\\X\5^% 

“ 1 “ 5 


aOaS; = a(')aS, aWaSI = -A(')aS, A<')aOT = f&A<')AS, 

^<t>y 

A'‘lAg = aI')a°t A®aJT = ^a(‘)aS, A(‘)aS = -A(‘>aS. 

^<f>y 

aWa^t =-A«AfT, A^A^^e =-A^A^T, A^^^ X^g = - A^H^. A^a^t = .aWa^^- 


D Computation of DR parameters at different scales 


The values of the DR parameters at the scale /U are obtained by renormalization group running 
from the starting scale hq to the scale /r. If the scales are not too far apart an approximate 
result can be obtained by exploiting the relation between OS and DR parameters p at the scale 

^os ^ Sp^^{p) = p^^ip) + Sp^^ip) . (D.61) 


Here 5p^^ and denote the OS and DR counterterm, respectively. The scale dependence 

in the DR counterterm, which purely subtracts the UV divergences, enters through the scale 
dependence of the parameters. As has been shown in 59 the only DR parameters that receive 
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two-loop counterterms at order 0{atas) are tan /3 and A, and this arises due to the non-vanishing 
wave function renormalization counterterm for Hu- We exemplify for tan/3 how to obtain the 
relation between the DR renormalized tan^S’s at two different scales /ri and /i 2 - We denote by 
tan/3P“®^^ the tan/3 defined through the DR condition with the top/stop sector renormal¬ 
ized DR. Analogously, tan/3P“®‘^® is understood to be the OS tan/3 and the top/stop sector 
renormalized in the OS scheme. The relation between these two definitions of tan /3 is given by, 

tan /3P'i'-«OS + ^(1) tan ^P^^eOS ^ ^(2) ^pureOS ^ 

tan^SP”*^^^ + (3(1) tan / 3 P“«DR ^ ^( 2 ) tan/3P”®°^ , 


where again the superscripts ( 1 ) and ( 2 ) refer to the one- and two-loop counterterm, respectively. 
The one- and two-loop counterterms in the pure OS scheme can be expanded in terms of e as 

.OS'! 

(D.63) 


(3(1) tan/3P“"OS ^^ 2 . 

(3(2) tan/3P""®°S = /" ( ^ a2(mps, oiDR(^)) 


+ /2(mP®,apii(^)) 1 (D.64) 


where the functions oi and fi do not depend on the renormalization scale ^ while 02 , 62 and /2 
implicitly depend on /x through their dependence on (^^^(/x). Note that these expansions can 
only be applied in the OS scheme of the top/stop sector in the context of our calculation. In 
the limit e —)• 0 these equations read 


(3(i)tan/3P“"°s = + ai(mP^) In/x^ +/i(mP^) (D.65) 

^( 2 ) ^pureOS ^ b2{m^^, a2(mP^apil(/x)) -F 2 62(mps,apil(/x))ln/x^ 


e 

os „,DR/ ,,2 I or X.v^OS ^,DR/..^M„2 ..2 


,OS „,DR 
L 

+2 02 (mP'’, a^ii(/x)) + 2 62 (ux^^, q;^ii(^)) In"' fi 

+/ 2 (mP®,apil(/x)) . 

The one- and two-loop counterterms in the pure DR scheme are 

,DIl 


e 

j(2)tan/3PureDR ^ b 2 {m^^{fi),a^^{fi)) ^ C 2 im^^{fi),a^^{fi)) 


(D. 66 ) 

(D.67) 

(D. 68 ) 


Replacing Eqs. ( |D.65[ ), ( |D. 66 [ ), ( |D.67[ ) and ( |D. 68 | ) into Eq. ( |D.62[ ) one gets the relation of the 
pure DR renormalized tan/3’s at the scales /xi and fi 2 - Taking into account the relation 


mP^ = mP® -b ((3mt)fin , (D.69) 

where ((3r?T,f)fin denotes the finite part of the OS counterterm, all terms proportional to the poles 
in e cancel at the considered order, and we are left with 


_ _ ,,2 

tan/3P“®°ii(/xi) - tan,5P“'^iiil(;U2) = ai(mpS)ln^ 

fJ -2 

+2 a2{mf^,a^^{fii)) In fil - 2 a2(mP®, Q;pil(/X2)) ln/X2 
+2b2{mf^, apil(/xi)) In^ fif - 2b2{m^^,a^^{fi2)) In^ fi^ - (D.70) 
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For the parameters p that are renormalized at one-loop order only, this relation simplifies to 

_ _ ,,2 

^pureDR(^^) _ ^pureDR^^^) = In « . (D.71) 

fJ'i 
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